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ABSTRACT 

Transiting circumbinary planets are more easily detected around short-period than 
long-period binaries, but none have yet been observed by Kepler orbiting binaries 
with periods shorter than seven days. In triple systems, secular Kozai-Lidov cycles 
and tidal friction (KLCTF) have been shown to reduce the inner orbital period from 
~ 10'^ to a few days. Indeed, the majority of short-period binaries are observed to 
possess a third stellar companion. Using secular evolution analysis and population 
synthesis, we show that KLCTF makes it unlikely for circumbinary transiting planets 
to exist around short-period binaries. We find the following outcomes. (1) Sufficiently 
massive planets in tight and/or coplanar orbits around the inner binary can quench the 
KL evolution because they induce precession in the inner binary. The KLCTF process 
does not take place, preventing the formation of a short-period binary. (2) Secular 
evolution is not quenched and it drives the planetary orbit into a high eccentricity, 
giving rise to an unstable configuration, in which the planet is most likely ejected from 
the system. (3) Secular evolution is not quenched but the planet survives the KLCTF 
evolution. Its orbit is likely to be much wider than the currently observed inner binary 
orbit, and is likely to be eccentric and inclined with respect to the inner binary. These 
outcomes lead to two main conclusions: (1) it is unlikely to find a massive planet on 
a tight and coplanar orbit around a short-period binary, and (2) the properties of 
circumbinary planets in short-period binaries are constrained by secular evolution. 

Key words: gravitation - celestial mechanics - planet-star interactions - stars: 
kinematics and dynamics. 


1 INTRODUCTION 

An increasing number of transiting circumbinary planets 
around solar-type main-sequence (MS) binary stars are be¬ 
ing discovered by the Kepler mission. Currently, 10 such 
planets are known (see Table \M\ top rows). The orbital 
periods of the stellar binaries in the systems discovered 
so far have a mea n value of 20.4 d; the shortest period is 
7.45 d (Kepler 47: IOrosz et al.l[2012al ). In contrast, the Ke¬ 
pler eclipsing binaries without circumbinary planets typically 
have shorter orbital periods, with a mean period of 2.8 d. 
There exists a bias for detecting more transiting circumbi¬ 
nary planets around shorter-period binaries because (1) the 
shorter period results in more transits in a given amount 
of time, and (2) the binary precession period is only a few 
years long, resulting in intersecting binary and planet orbits 
(as seen in the plane of the sky) during the Kepler mis- 
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sion dMartin &: Triau3 l2015l ). Because of this bias, many 
circumbinary planets are expected to have been observed 
around short-period Kepler eclipsing binaries. However, so 
far, none have been found. 

If the apparent lack of (nearly) coplanar circumbinary 
planets around close binaries is intrinsic and not related to 
(yet unknown) observational bias(es), then this raises the 
question of its origin. The following two observations suggest 
that this origin is related to the (MS) evolution of solar-type 
triple star systems. 

• The fraction of tertiary companions to spectroscopic 
binaries is a strong function of the inner period, increasing 
from 0.34 for inner periods Pa > 12d to 0.96 for Pa < 3d 
llTokovinin et al.ll2006h . For the Kepler eclipsing binaries in 
particular, the tertiar y fraction for inner binaries with pe- 
riods < 3d is ~ 20% (iRappaport et al.ll2013l : IConrov et al.l 
I2OI4I) . The latter studies are limited (by observing time) 
to triples with outer periods < 3yr. The complete tertiary 
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fraction for these Kepler binaries, i.e. including any outer 
period, is likely much larger. _ 

• In the survey of iTokovininI ll2014h . a peak is found in 
the (inner) period distribut i on around 3 d. A s show n by 
iFabrvckv fc Tremaine! (l2007f ): lNaoz fc Fabrvckvl ll2014h . this 
peak in the period distribution can be explained by the com¬ 
bination of the secular gravitational torque of the tertiary 
companion and tidal friction in the inner binary . The f ormer, 
throu gh Kozai-Lidov (KL) oscillations dLidovI 1 19621 : iKozail 
Il962h . can excite the inner binary orbit to high eccentric¬ 
ity. Consequently, tides are much more effective, leading to 
tidal dissipation and shrinkage of the orbit, and, eventually, 
to a tight and (nearly) circular orbit. In this scenario, the 
precursors of the binaries affected by KL cycles with tides 
have orbital periods ranging roughly between 10 and 10^ d. 

The lower limit of 10 d of precursor binaries is right in 
the ballpark of the orbital periods of the currently observed 
Kepler systems with transiting circumbinary planets. Also, 
there are a number of observed wider systems orbited by at 
least one circumbinary planet; some of them are listed in 
Table IXT] (bottom rows). Therefore, a large fraction of these 
precursor binaries might have a circumbinary planet in a 
stable orbit in-between the ‘inner’ and ‘outer’ stellar orbits 
(henceforth, we refer to the inner and outer stellar binaries 
simply as the ‘inner’ and ‘outer’ binaries; see also Fig. [T]). 

Here, we show that in these stellar triples with a cir¬ 
cumbinary planet around the inner binary, the planet can 
strongly affect the secular orbital evolution of the inner bi¬ 
nary compared to the situation without a planet, provided 
that certain conditions are met (see below). This ‘shielding’ 
effect arises from a quenching of the eccentricity oscillations 
in the inner binary induced by the torque of the outer bi¬ 
nary, because of precession induced from the circumbinary 
planet. The shielding effect is similar to the quenching of KL 
oscillations in stellar triples because of additional sources of 
orbital precession (notably, precession associated with gen¬ 
eral relativity, tidal bulges and/or stellar rotation). In some 
cases, however, resonances can occur that enhance the ec¬ 
centricity oscillations in the inner binary (such resonances 
can also o ccur in isolated trip les in conjunction with general 
relativity; iNaoz et al.ll2013'blf . 

Typically, shielding of KL oscillations by the presence 
of the planet is only effective for massive planets that are 
in a (nearly) coplanar and tight orbit with respect to the 
inner binary. This implies the following two scenarios for 
the secular evolution of the system. 

• If the precursor binary is orbited by a massive planet in 
a tight and coplanar orbit, then planet shielding can prevent 
the binary from shrinking through the combined effects of 
KL cycles and tidal friction (KLCTF). 

• On the other hand, if the planet is of low mass and in 
an inclined and wide orbit with respect to the inner binary, 
then shielding is typically ineffective, and the binary can 
shrink through KLCTF to become a short-period binary. 

Assuming that short-period binaries are produced through 
KLCTF, then this implies an intrinsic lack of massive plan¬ 
ets in tight and coplanar circumbinary orbits (i.e. strongly 
shielding planets), whereas far-away and inclined low-mass 
circumbinary planets (i.e. weakly shielding planets) could 
still be abundant. This is consistent with the current null¬ 
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Figure 1. A schematic depiction (not to scale) of the hierarchi¬ 
cal configuration of the systems considered in this paper. The 
circumbinary planet (m 3 ) orbits the inner binary (mi and m 2 ); 
the outermost star (m 3 ), in the outer binary, orbits the centre of 
mass of the inner binary-fplanet system. 


detections of transiting Kepler circumbinary planets around 
short-period binaries. 

The goal of this paper is to quantify the above argu¬ 
ment. The organization is as follows. In Section[2l we briefly 
describe our methods and assumptions. In Section|3l we 
study the secular gravitational dynamics of stellar triples 
with a circumbinary planet, and we quantify the conditions 
when the planet can affect the inner binary. In order to 
evaluate the effect of circumbinary planets in a population 
of triples in the field, we carry out a population synthesis 
study in Section[4l We discuss our results in Section[5l where 
we also describe, using direct A-body integrations, the fate 
of circumbinary planets that become unstable because of 
secular evolution. We conclude in Section[6] 

We remark that nearing the completion of this pa¬ 
per, we became aware of two similar and independent 
studies on the lack of circumbina ry planets around s hort- 
period binaries in st ellar triples bv iMartin et al.l (l2015l f and 
iMunoz fc Lail ll2015l '). 


2 METHODS AND ASSUMPTIONS 

We model the system of a circumbinary planet orbiting 
the inner binary in a stellar triple system as a hierarchi¬ 
cal quadruple system in the ‘triple-single’ c onfigu ration, 
which we studied previously in iHamers et al.l ll2015h (h ere- 
after HPAPZ15). In the latter work, the Hamiltonian of the 
system was derived and expanded in terms of the ratios of 
the three binary separations rA, re and rc, where, by as¬ 
sumption, rc re rA. In the current context, binary A 
corresponds to the inner binary, binary B to the orbit of the 
planet around the inner binary, and binary C to the orbit 
of the tertiary star around the centre of mass of the inner 
binary-bplanet system (evidently, the latter nearly coincides 
with the centre of mass of the inner binary). For consistency 
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with HPAPZ15, we use indices 1 and 2 to denote quantities 
associated with the inner binary primary and secondary, re¬ 
spectively, index 3 for the planet, and index 4 for the tertiary 
star. 

A schematic depiction of our configuration is shown in 

Fig-IIl 

In HPAPZ15, the orbit-averaged Hamiltonian was de¬ 
rived, and a numerical algo rithm was developed within 
the AMUSE framework dPorteeies Zwart et al.l l2013l : 
IPelupessv et al.l l2013h to solve the resulting equations of 
motion. Post-Newtonian (PN) dynamics at the IPN or¬ 
der are taken into account as described by equation (9) of 
HPAPZ15. 

Here, we extended this algorithm by also including the 
effects of tidal friction within the inner binary. Gravitational 
perturbations by the planet and the tertiary star with re¬ 
gards to the Inner binary tidal evolu tion were ignored. W e 
adopted the equilibrium tide model jEggleton et al.lll998ll . 
in which it is assumed that each star quasi-hydrostatically 
adjusts its shape to the time -varying perturbing pote ntial 
of its companion. Following iBarker fc Ogilvia (l2009l l. we 
adopted a constant tidal quality factor Q related to the 
mean motion n and the tidal lag time r via Q — l/(nr). 
Below, instead o f Q, we use the directly related quantity 
Q' = 3Q/{2k) llBarker fc Qgilvij l2009l l. where k is the 
second-order potential Love number. Typical values of Q' 
for solar-type stars , as inferred from observations, are Q' ~ 
5.5 X 10® — 2 X 10® IIMeibom &: Mathieijl2005| : lOgilvie fc LirJ 
bOOit) . 

A constant Q implies that r = l/(nQ) = P/{2irQ) 
scales with the orbital period, and, therefore, r effectively 
decreases as the orbital period decreases due to tidal fric- 
tion. This may not give an e ntirely accurate descr iption; e.g. 
ISocrates fc Kat j (l2012f ) and lSocrates et ^ ll2012ll show that 
with a numb er of simplifying assumptions, r is constant, and 
the results of lHanserJ ll2010ll suggest that Q' is not constant 
and depende nt on the semimajor axis ( cf. equation 13 of 
lHaJiserJ l2010l l. Nevertheless, as argued by [Barker fc Qgilvid 
(l2009tl . given the current uncertainties in the underlying 
tidal dissipation mechanisms and, therefore, the efficiency 
of tidal dissipation, the assumption of a constant Q is useful 
for studying the general effects of tidal friction. 

The equilibrium tide model is included in our algo¬ 
rithm by adding the relevant terms deA/dt|TF, dhA/dt|TF, 
df2i/dt|TF and dJ 22 /dt|TF in th e equations of mo t ion, a s 
given by equations (A7)-(A15) of I Barker fc Qgil-^ ll2009ll . 
Here, ej, and hk are the eccentricity vector and the specific 
angular momentum vector of orbit k, respectively, and flk 
is the spin angular momentum vector of star k £ {1,2}. 
Throughout this paper, the initial spin vectors flk are as¬ 
sumed to be parallel with the inner orbital angular momen¬ 
tum vector fiA- 

Our algorithm, that models the orbital evolution, 
is coupled within AMUSE with the ste l lar evolution 
code SeBa IPortegies Zwart fc Verbuntlll99^ : lToonen et al.l 
I 2 OI 2 I I. which is also interfaced within AMUSE. We use the 
latter code to compute the masses and radii during in the 
integration. In this work, we focus on low-mass MS stars. 
Therefore, the time-dependence of the latter quantities is 
weak for time-scales less than a Hubble time. We remark, 
however, that this is no longer the case for post-MS evolu¬ 
tion. 
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Figure 4. Illustration of the joint dependence of the maximum 
eccentricities in binaries A and B on a a and 13 , as an extension 
of the top-middle panel of Fig. [2] In each panel, 1 — e^ax k 1® 
plotted as a function of i^. The value of ae is indicated for each 
panel; for all panels, m 3 = 1 Mj. The black solid horizontal lines 
show the value of 1 — Gmax.A in the absence of the planet. 


3 THE PLANET-SHIELDING EFFECT 

Before studying the effect of circumbinary planets in a pop¬ 
ulation of triples (cf. Sectlon|4]), we hrst focus on a slmplihed 
case where the stellar triple system is kept fixed, and plan¬ 
ets are ‘Inserted’ at stable orbits at various loci in-between 
the inner and outer binaries. The aim is to gain quanti¬ 
tative insight into the effect of the planetary orbit on the 
eccentricity oscillations and the tidal evolution In the inner 
binary. For simplicity, we here focus on triples in which the 
octupole-order terms are unimportant. Note, however, that 
the (non-cross) octupole-order terms are always included in 
the integrations (cf. Section [531l . and no restriction is made 
with respect to the importance on the octupole-order terms 
in Section|4l 


3.1 Kepler transiting circnmbinary systems 

Here, we focus on systems similar to the Kepler transit¬ 
ing circumbinary systems. The inner semimajor axes in the 
latter are a a ~ 0.1 — 0.2 AU (cf. the top rows of Table 
mj, and relativistic precession is therefore important in 
these systems . As is well known from previous studies (e.g ., 
iHolman et al.l 1 19971 : iBlaes et"^ I 2 OO 2 I : iNaoz et ^ 12013^ , 
general relativistic precession limits the range of semima¬ 
jor axes and eccentricities of tertiary orbits for which the 
inner orbit eccentricity is excited. In general, an order-of- 
magnitude estimate of the KL time-scale for the binary pair 
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Section_ O_O_O_S 


mi/M0 

1 

1 

0.83 

0.5-1.2 

m2/M0 

0.5 

0.5 

0.50 

0.1-1.2 b 

m^lMj 

0.01-1 

0.001-3.16 

0.01-1 

0.01-1 

rtii/MQ 

1 

1 

0.73 

0.1-1.2 

oa/AU 

0.2 

1 

0.20 

0.03-2.0 

ob/AU 

1-5 

10-50 

1.09 - 1.22 

0.14- 





3.1x102 ® 

ac/AU 

20 

160 

7.73 

1.5- 





1.8x10-* 

ea 

0.01 

0.01 

0.40 

0.01-0.95 f 

SB 

0.01 

0.01 

0.01 

0.01 

ec 

0.01 

0.01 

0.30 

0.01-0.95 f 

iAP 

0.01 

0.01 

0.01 

0.01 

iBp 

0.01-180 

0.01-180 

0.01-180 

0.01-180 

ic/° 

75 

85 

73.7 

0.01-180 s 


0.01 

0.01 

287.7 

0.01-180 

‘*Jb/° 

0.01 

0.01 

281.0 

0.01-180 ** 

‘^c/° 

0.01 

0.01 

51.6 

0.01-180 

UaP 

0.01 

0.01 

166.1 

0.01-180 

flBp 

0.01 

0.01 

42.6 

0.01-180 

ncP 

0.01 

0.01 

340.1 

0.01-180 ** 

^spin,fe/d 

N/A 

10 

10 

10 

Q'k/W^ 

N/A 

1 

1.9 

0.55-2 h 


N/A 

0.08 

0.08 

0.08 

Salpeter distribution dN/dmi oc 

Sampled from m2 = where 0 < gin < 1 is linearly distributed. 

Sampled from m^ = gout (mi -|- m2), where 0 < gout < 1 is linearly 


distributed. 

^ Lognormal distribution in the orbital period with mean 

logio(Pfc/d) = 5.03, standard deviation = 2.28 and range 

—2 < log]^o(^fc/d) < 10. Triples are subject to stability constraints, and 
the inner binary should not merge in the absence of a planet. 

® Linearly sampled for each triple, with lower limit i = 1.5 UB,crit,ABi 
where crit AB is the critical semimajor axis for dynamical stability 
of the planet in a coplanar orbit around the inner binary, and which i s 
adopted from the fitting formula given by llHolman Sz WiegertlllOOil^ . 
The upper limit is a.B,u = 0.9 ub, crit,BC > where aB,crit,BC i® fh® largest 
possible value of ub lor dynamical stability with respect to the orbit of 
the tertiary star, estimated by applying the lMardling Aarsetl]l20Qll 
criterion to the BC pair, with an ‘outer’ mass ratio of gout = m^Kmi + 
m 2 + m3). 

^Sampled from a Rayleigh distribution, dN/de^ oc exp(—/3e|), with 

rms = 0.33 llRaghavan et al.ll2010l) . 

s Sampled from a linear distribution in cos(2c), with —1 < cos(2c) < 1- 
^ Sampled from a linear distribution. 

Table 1. Initial conditions for the systems in Section r3 .1 1 (first column, 
cf. Figs [2l4t . for the systems in Section[3]2] (second column, cf. Figs [5] 
an dEll, for the systems in Section FS.31 fthird column, cf. Fig. l20|l . and for 
the population synthesis in Section[4| (fourth column). Note that for the 
third column, the systems are integrated with the secular code until the 
orbit of the planet intersects with the inner binary; subsequent evolution 
is carried out using direct N-body integration (cf. Section r5.3ll . The or¬ 
bital angles (inclinations arguments of pericentre cok and longitudes 
of the ascending node are defined with respect to an arbitrary, but 
fixed reference frame; see equation (6) of HPAPZ15 for the relation be¬ 
tween these angles and the orbital vectors. The spin periods Pspin,fei 
tidal quality factors and radii of gyration rg^fc apply to the inner 
binary stars and are used for tidal evolution. 


kl is given by 


Pklm 


Pi mk,p + mk,s + mi,'. 

Pk 


(1 


2 \ 3/2 

-ei) 


where mk,p = mi, mk,a = m 2 and mj^g = m 3 in the case of 
-Pkl.ab, mk,p = mi + m2, mk,a = m3 and mi,s = m^ in the 
case of Pkl.bc, and m^.p = mi, mk a = m 2 and mi^a = mi in 
the case of Pkl. ac (linnanen et al. I ll997l . HPAPZ15, see also 
lAntogninll2015l h The (pairwise) IPN time-scale for orbit k 


© 2015 RAS, MNRAS 000.ITH241 















(^Ab)/° ^ 6max,A: 


Circumbinary planets in triples 5 


Qo 

0.13 1.00 7.63 58.24 



Figure 2. Top row: maximum eccentricities in the inner binary (red bullets) and circumbinary orbit (blue stars) for the integrations of 
Section FS-ll fcf. the left column of Table [J . as a function of the initial as (left column), (middle column) and m3 (right column). In 
the left column zb = 0° and m3 = 1 Mj; in the middle column = 1-2 AU and m3 = 1 Mj; and in the right column aB = 1.2 AU and 
Zb = 0°. The corresponding value of Qq (cf. equation [8ll is indicated in the top left panel. In the bottom row, the average inclination 
between the inner and circumbinary orbits is shown. The two black vertical dashed lines in the top right panel indicate Earth and Jupiter 
masses, respectively. 


is given by 

tiPN,k = -Pk (1 — el) —( 2 ) 
3 rg,fe 

where rg^k = Gmtot,fc/c^, with mtot.A = mi + m 2 , mtot.B = 
mi + m 2 + m 3 and mtot,c = mi + m 2 + m 3 + m 4 , is the 
gravitational radius. Equating equation m as applied to the 
AC pair, to equation as applied to binary A, we find that 
relativistic precession dominates in binary A if 

(1 — eX) m 4 (mi + m 2 + m 3 + m 4 ) 

3G(mi + m 2 )^(mi + m 2 + 1714 ) 
SS34AU, (3) 


ac > (1 - ec) 


An estimate of the lower limit on ac for dynami¬ 
cal stability with respect to the inner binary (not yet 
taking into account the circumbinary planet, and ignor- 
ing the dependence on m utual inclination) is given by 
llMardling fc Aarsethll200ih 

2.8aA 


(4) 


1 -f 


m4 


1 -f ec 


mi -I- m2/ VI - ec. 


: 0.67 AU. 


1/3 Here, we assume a fixed ac = 20 AU. With this choice, the 
octupole parameter toct.AC, defined as 


eoct.AC 


mi — m2 oa ec 
mi -I- m2 ac 1 — Cq’ 


(5) 


where for the numerical estimate we assumed ca = ec = 
0.01, mi = 1 Mq, m 2 = 0.5 M©, m 3 = 1 Mj and mj = 1 Mq. 
In the remainder of Section ISTTI we shall assume these values 
for the numerical estimates. 


is Coct.AC ~ 1.7 X 10“®, indicating that octupole-order terms 
in the Hamiltonian expansion for the AC pair are unimpor¬ 
tant |Li thwick &: Naod I 2 OIII: IXatz et al]l 201 ll : iNaoz et al.l 
l2013ai : Tevssandier et al.] 20131 : Li et al.l 20l4 l. Assuming an 
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Q Q Lk-T .> _^_ . . . t ..i. L-<r? ' . I. . ■ >. • . . I ■ ■ . ! , I 

0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 


t/Myr t/Myr t/Myr 

Figure 3. Illustration of the time evolution of the eccentricities of the inner orbit and the circumbinary, for six systems corresponding 
to the top-left panel of Fig. (2] In each panel, the value of Qo is indicated. 


initial mutual inclination of iac = 75° between binaries A 
and C, and with the inclusion of the IPN terms, the maxi¬ 
mum eccentricity attained in binary A, without a circumbi¬ 
nary planet, is CA.max ~ 0 . 88 . 

As a next step, we include a circumbinary planet with 
mass ms in orbit around binary A with semimajor axis 
ttB < ac, and with various mutual inclinations iab. For 
simplicity, the circumbinary orbit is assumed to be initially 
nearly circular, i.e. es = 0.01. Assuming iab = 0°, for dy¬ 
namical stability with respe ct to the inner binary, the semi¬ 
major axis Ob must satisfy dHolman fc Wiegerd 19991 ) 

flB > aA/(eA,MA) ~ 0.48 AU, ( 6 ) 


where /(ca, /ia) is a function of the inner binary eccentricity 
ca an d mass ratio ha = m 2 /(m i -|- m 2 ), given by equation 
(3) of lHolman fc Wieeertl (1 19991 ). An estimate of the largest 
possible value of as for dynamical stability with respect to 
the outer orbit, bi nary C, ignoring dependen cies on inclina¬ 
tions, is given by jMardling &i Aarsethll200 j ) 


Ob < ac 


1 — ec 


2.8 

5.8 AU. 


1 + 


m4 


mi -h m 2 -I- m 3 y ^1 — ec. 


1 + ec 


1 -2/5 


(7) 


Note that, strict ly _ speaking, the criterion of 

llMardling fc AarsethI I 2 OO j) only applies to three-body 
systems. However, by carryi ng out a number of A- body 
integratio ns with Mikkola llMik kola fc MerrittI l2008l) and 
Sakura llGoncalves Ferrari et al.l 20141 ) within AMUSE, 
we find that circular circumbinary orbits with ab = 5.8 AU 
in the system of our choice are indeed dynamically stable 
on time-scales of a few multiple periods of the outer binary. 

Based on these estimates, we carried out integrations of 


the long-term secular dynamical evolution with our orbit- 
averaged algorithm (cf. Section[21), with ob ranging between 
1 and 5 AU, iab between 0 and 180°, and m 3 between 10~^ 
and 10° Mj. To simplify the interpretation, tidal evolution 
is not included in Section lXB The integration time was set 
to 2Pkl,ac, i.e. the time-scale for two KL oscillations in the 
AC pair, without taking into account the effect of the planet. 
A comprehensive list of initial parameters is given in the left 
column of Table [T] 

In the left, middle and right panels of the first row of 
Fig. [ 2 ] we show the resulting maximum eccentricities in bi¬ 
naries A (red bullets) and B (blue stars) as a function of 
Ob (with iB ~ 0° and m 3 = IMj), ib (with ae = 1.2 AU 
and m 3 = 1 Mj), and m 3 (with ob = 1.2 AU and zb = 0°), 
respectively. In the bottom row of the same figure, we show 
the average inclinations between orbits AB as a function of 
Ob, *b and m 3 . In the top left panel, we also show an equiva¬ 
lent dependence on the parameter Qo, defined as the initial 
ratio of the KL time-scales for the AB and AC pairs. 


_ Pkl.ab.o 

Qo = - 

Pkl.ac.o 


flB 

ac 


I , , I 

m 4 mi -I- m 2 -I- m 3 -I- m. 
m 3 mi -I- m 2 -I- m 4 


/rru f 1 - eB,o \ 
\acj rnsyi-eloj 




zii 


(8) 


Here, in the third line, we neglected the planet mass com¬ 
pared to the sum of all stellar masses. We note that Qo is 
closely related to the quantity Po defined by HPAPZ15, and 


© 2015 RAS, MNRAS 000.[TIf24l 
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is given by 


TZo = 


fkL.AB.O 

ft:L,BC,0 

VoAfl^y 

\a.AalJ 


Qo(^ 

aA 


mi + m2 
mi + m 2 + m 3 


1/2 


m4 

m 3 


1 ~ ^1,0 
1 “ ®c,o 


m4 

m 3 


1 ~ ^B ,0 
1 “ ®C,0 


3/2 


3/2 


3/2 


(9) 


If iB = 0° and Qo ^ 1, fimax.A ~ 0.88 is the same 
as in the case without the planet (cf. the top-left panel of 
Fig. m- In this limit, the torque of the outer orbit on the 
inner orbit (for which is a proxy) is much larger 

than the torque of the circumbinary orbit on the inner orbit 
(for which Pj^^ ab ^ ac i® P^'O^y)- The circumbinary 
orbit therefore only has a small perturbative effect on the 
inner orbit, and the maximum eccentricity is not affected. 
If Qo 1, then also TZo S> 1 (cf. equation [9|. This implies 
that the torque of binary B is not strong enough to maintain 
coplanarity between binaries A and B, as is shown in the 
bottom-left panel of Fig. (2] Effectively, binary B decouples 
from binary A, and the torque of the outer orbit results in 
high-amplitude eccentricity oscillations in binary B. 

The latter oscillations have large enough eccentricities 
for the circumbinary orbit to intersect with the inner bi¬ 
nary, i.e. rp,B = nB(l —cb) < fa ,a = aA(l + eA). This is the 
case for all points below the black solid line in the top-left 
panel of Fig. [21 For the purposes of this section, no stop¬ 
ping conditions were imposed on the eccentricities during 
the integrations. However, in reality, the circumbinary orbit 
is dynamically unstable for points below the black solid line. 
This is borne out by integrating some of the systems below 
this line with direct A-body integration, where the orbital 
parameters are set to correspond to the moment of maxi¬ 
mum eccentricity of binary B. We investigate the possible 
outcomes for these cases of dynamical instability further in 
Section lAdl 

On the other hand, if is = 0° and Qo ^ 1, emax,A ~ 
eA,o. In this case, the net torque on the inner orbit is dom¬ 
inated by the circumbinary orbit. Consequently, the inner 
orbit processes much more rapidly compared to the case 
without the planet, and this results in a quenching of the 
KL eccentricity oscillations otherwise induced by the outer 
orbit. Furthermore, the period of the latter oscillations is 
substantially reduced, as is illustrated in Fig. O where we 
show the time evolution of ba and bb for several values of Qo, 
assuming ib = 0° and m 3 = 1 Mj. Note that the inclination 
between orbits A and B remains zero in this case (cf. the 
bottom-left panel of Fig. o, and that octupole-order terms 
are not dominant, even for our smallest value ab = lAU 
(for flB = lAU, Eoct.AB ~ 6.7 X Therefore, the cir¬ 

cumbinary orbit, on its own, does not induce eccentricity 
oscillations in binary A. 

Whenever the circumbinary orbit becomes unstable, the 
planet could collide with one of the stars, or be ejected from 
the system (see also Section lS^ . Evidently, in either case, 
the planet can no longer shield the inner binary from KL- 
eccentricity oscillations induced by the outer orbit. Typi¬ 
cally, this occurs whenever the circumbinary orbit is rela¬ 


tively close to the outer orbit (Qo > 1), and is therefore 
strongly affected by the latter orbit’s torque. In the same 
regime, however, shielding is ineffective. This implies that 
dynamical instability does not necessarily rule out the pos¬ 
sibility for shielding by the planet. In the top-left panel of 
Fig. 121 this is indeed the case: the circumbinary orbits be¬ 
come unstable for Qo ^ 1, whereas shielding only occurs for 
Qo < 1. 

The dependence of the maximum eccentricities on ib is 
shown in the top middle panel of Fig. [2] for ub = 1.2 AU. 
There is a strong dependence of the planet’s shielding ability 
on iB. For coplanar (either prograde or retrograde) orbits, 
eA.max is closB to the initial value ca.o- As the initial orbit of 
the planet becomes more inclined with respect to the inner 
binary, shielding becomes less efficient. Interestingly, for ib 
close to 90°, shielding is again more efficient. 

The nature of the dependence on ib also strongly de¬ 
pends on flB. In Fig. 0] we show the dependence on ib for 
other values of ob, assuming m 3 = IMj. Generally, Cmax.A 
is weakly dependent on ib for either small ab or Qq ^ 1 
(in which case emax.A ~ ca.o), or large ob or Qo » 1 (in 
which case Cmax.A is approximately equal to the value in ab¬ 
sence of the planet). In the intermediate regime, there is a 
complicated dependence on Tb- Generally, quenching of the 
KL eccentricity oscillations in binary A is most efficient for 
close to coplanar orbits (in the case of Fig. (4] 0° < ib 30° 
and 150° < *b < 180°). Furthermore, for some values of 
Ob, quenching can also be effective for ib ~ 90° (cf. the top 
middle panel of Fig. [2] and the top left panel of Fig. 0]). 

The dependence on m 3 is illustrated in the top-right 
panel of Fig. [2] (for ob = 1.2 AU and zb = 0°). The planet 
mass must be large enough for shielding to be effective, 
which is intuitively easy to understand. In our example, a 
Jupiter-mass planet can effectively shield the inner binary, 
assuming that it is coplanar with and close to the inner bi¬ 
nary, whereas an Earth-mass planet cannot provide shield¬ 
ing, regardless of its orbit. 


3.2 A triple with a shrinking inner binary orbit 

Having considered in Section ISTi the planet’s ability to shield 
the inner binary from KL eccentricity oscillations induced by 
the outer orbit, we here extend the analysis by also including 
the effects of tidal friction in the inner binary. We choose a 
triple system with a relatively wide inner binary (oa = 1 AU 
or Pa ~ 300d) and tight outer binary (ac = 160 AU), such 
that in the absence of a planet, the inner binary becomes 
highly eccentric because of KL eccentricity cycles induced 
by the tertiary, and shrinks to a tight binary with Pa ~ 8 d. 

As in Section [3Tl we carry out a set of integrations with 
a planet with various values of ob, *b and m 3 . A compre¬ 
hensive list of the initial conditions is given in the second 
column of Table [T] For each system, the integration time is 
set to IOPkl.ac ~ 50Myr. In contrast to Section01 here we 
stop the integration if the circumbinary orbit intersects with 
the inner binary, i.e. if rp,B = aB(l — cb) < ra,A = aA(l-l-eA) 
(neglecting the radii of the stars and the planet). We assume 
that in the latter case, the planet either collides with one of 
the stars or is ejected from the system (see also Section iTSll . 
and, therefore, the subsequent evolution is equivalent to the 
situation without a planet, i.e. the inner binary shrinking to 
an inner final inner period Pa,! « 8 d. 
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Figure 6. The final inner orbital period for the integrations of Section|3]2] (cf. the second column of Table[TJ as a function of the initial ob 
( left column), ie (middle column) and m 3 (right column). In the left column jb = 0° and m 3 = 1 Mj, in the middle column ob = 40 AU 
and m 3 = IMj, and in the right column ub = 40 AU and ib =0°. Cases where the planet became unstable are denoted with black 
crosses. 


In Fig. [3 we show, as an example, the evolution of the 
inner binary semimajor axis a a (top panel), the inner bi¬ 
nary eccentricity ba (bottom panel; black lines), and the cir- 
cumbinary eccentricity cb (bottom panel; red lines), where 
Ob = 40 AU and ib = 0° initially. Two planet masses are as¬ 
sumed: m 3 = 10“® Mj (solid lines) and m 3 « 1.32 x 10“^ Mj 
(dashed lines). The values of Qo are ~ 55 and « 0.42 for 
these values of m 3 , respectively. 

For the low planet mass (Qo « 55), shielding is ineffec¬ 
tive, and the inner orbit eccentricity becomes highly excited, 
resulting in efficient shrinkage of the orbit. Interestingly, this 
also increases the relative strength of the torque of the outer 
orbit compared to the inner orbit on the circumbinary or¬ 
bit. Consequently, the eccentricity of the circumbinary or¬ 
bit becomes excited at t ~ 25 Myr. Although not the case 
here, the latter eccentricity could become high enough for 
the circumbinary orbit to intersect with the inner binary, 
and become unstable. In other words, KL cycles with tidal 
friction tend to destabilize the circumbinary planet, and this 
further reduces the chances of a circumbinary planet resid¬ 
ing in a stable orbit around a binary that shrank due to KL 
cycles with tidal friction. This consequence of inner binary 
shrinkage on the pl anet orbital stability was also noted by 
iMartin et al.l (l2015l ) and iMunoz &: ll2015l ). 

In contrast, for the high planet mass (Qo « 0.42), 
shielding is effective, and the inner binary eccentricity does 
not become high enough for efficient tidal friction. Never¬ 
theless, ca still oscillates with an amplitude of « 0.6. The 
circumbinary orbit eccentricity does not become excited in 
this case. 

In the left, middle and right panels of Fig. ( 6 ] we show 


the final inner period PA,f as a function of ub (with ib = 0 ° 
and m 3 = lAfj), is (with as = 40 AU and m 3 = IMj), 
and m 3 (with ab = 40 AU and is =0°), respectively. Cases 
where the planet became unstable are denoted with black 
crosses. There is a strong dependence of the final inner pe¬ 
riod on the circumbinary parameters. If ub is sufficiently 
small (the regime Qo <C 1), shielding is effective, and the 
inner binary does not shrink. For larger ub, the planet be¬ 
comes unstable, and is therefore unable to shield the inner 
binary. Note that in this regime of Qo > 1, shielding would 
have been ineffective even if the circumbinary orbit would 
be stable (cf. Fig. [2|. 

Similarly, there is a strong dependence on is ■ For copla- 
nar (either prograde or retrograde) circumbinary orbits, 
shielding can be effective (if ab is small enough) and the 
inner binary does not shrink. However, for highly inclined 
orbits, the planetary orbit is unstable and, subsequently, the 
inner binary shrinks. The boundaries between these regimes, 
iB ~ 45° and ib ~ 135°, are very similar to the critical val¬ 
ues for eccentricity excitation due to the KL mechanisms 
(assuming quadrupole order). We note that, for much more 
compact circumbinary planets with a small ratio aB/oA, 
instability can occur for even coplanar systems lILi et al.l 
l2014h . 

For the value of ub and ib in the right panel of Fig. 
[ 6 l the planetary orbit is stable for all assumed masses m 3 . 
Effective shielding, and hence shrinkage of the inner orbit, 
can only occur if the planet is sufficiently massive, in this 
case, if m 3 > 0.1 Mj. 
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Figure 5. Example evolution of the inner binary semimajor 
axis (top panel), the inner binary eccentricity (bottom panel; 
black lines), and the circumbinary eccentricity (bottom panel; 
red lines), for the triple system discussed in Section|3]2] Two 
planet masses are assumed: m 3 = 10 ~® Mj (solid lines) and 
m 3 Ri 1.32 X 10“^ Mj (dashed lines). In either cases, 03 = 40 AU 
and jb = 0 ° initially. 


4 PLANETS IN TRIPLES WITH 

SHORT-PERIOD INNER BINARIES 

4.1 Initial conditions 

In Section[3l we demonstrated the planet’s ability to shield 
the inner binary assuming a fixed triple system, and used 
specific worked-out examples to understand the evolution of 
circumbinary planets affected by a third companion. In the 
following, we study the effect of circumbinary planets in the 
more general case, and consider a large orbital phase-space 
through the study of a population of triples, for which the 
inner binary, in absence of a planet, would be affected by 
KL cycles with tidal friction. 

For s tellar triples, a number of surveys of solar-type 
MS stars ( Duouem^ fc_^dayoj[l99ll : Tokovinin et ahlBood : 
iRaghavan et al.l 2010l : iTokovinInl 20141 ) have given insight 
into the initial distributions of the masses and orbital pa¬ 
rameters. However, initial distributions for a circumbinary 
planet in orbit of the inner binary, for which observations 
are currently strongly biased, are very poorly constrained. 
Therefore, we generate initial conditions using a combina¬ 
tion of Monte Carlo sampling for the stellar triple system, 
and fixed grids for the circumbinary planet. Our procedure 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 


logio(^A/d) 


Figure 7. The initial (dotted lines) and final (solid lines) in¬ 
ner period distributions for the sampled triples as described in 
Secti 0 n r 4 .ll Top row (black lines): isolated binary evolution, i.e. 
taking into account only tidal friction. Middle row (blue lines): 
isolated triple evolution, i.e. also taking into account the torque 
of the outer orbit. Cumulative distributions for both cases are 
shown in the bottom panel. 


consists of the following four steps (a succinct summary is 
given in the fourth column of Table [T]). 

1. A stellar trip l e sys t em is sampled, s i milar ly to 
iFabrvckv fc Tremaind (l2007l ) ; iNaoz fc Fabrvck^ (l2014l ) . For 
the inner binary, the primary mass mi is sampled be¬ 
tween 0.5 and 1.2 M;?) assuming a Salpeter distribution, 
dN/dmi oc llSalpeterl 1 19551) . The secondary mass 

m 2 is sampled from m 2 = Qinmi with the constraint 0.1 < 
m 2 /M 0 < 1 . 2 , where gin is sampled from a fiat distribution 
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loglO (/shrink) 


Figure 8. The factor /shrink = o,A,i/o,A,l with which the in¬ 
ner binary semimajor axis shrinks for the sampled triples as de¬ 
scribed in Secti 0 n r 4 .ll Black (blue) lines correspond to isolated 
binary (triple) evolution. The black vertical dotted line indicates 
the value of /shrink chosen for the selection of systems with sub¬ 
stantially shrinking inner binaries. 


with 0 < gin < 1. The mass of the tertiary star, 7714 , with 
0.1 < m 4 /M 0 < 1 . 2 , is sampled from 7714 = {mi + m 2 )gout, 
where Qout is sampled from a flat distribution 0 < gout < 1 - 
The stellar radii, Ri, R 2 and Ri, are computed initially, 
and at subsequent times, using the stellar evolution code 
SeBa as described in Section^ The spins and fl 2 of the 
stars in the inner binary are initially assumed to be paral¬ 
lel with the inner orbital angular momentum, and the spin 
periods are assumed to be 27r/||$lfc|l = lOd. The tidal qual¬ 
ity factors Qfc for the inner binary primary and secondary 
stars (cf. Section[5]) are sam pled linearly between 5.5 x 10® 
and 2 x 10®, adopted from lOgilvie fc LinI l2007l . who con- 
strained Qj using observat ions of spectroscopic binaries by 
iMeibom fc MathieuI (l2005l ). The gyration radii for the 

primary and secondary stars are assumed to be = 0.08. 

The inner and outer orbital periods. Pa and Pc, are 
both sampled from a lognormal distribution with mean 
logj^Q(Pfc/d) = 5.03, standard deviation n'iogiQ(Pfc/d) = 2.28 
and range —2 < logj^g(Ps,/d) < 10. The latter distri¬ 
bution is co nsistent with the orb i tal p er iods of binaries 
as fou n d by Duquennoy fc Mavoil (Il99ll ): iRaghavan et al.l 
(l2010l) : iTokovininI ||2014|) . We choose to impose the fur¬ 
ther restriction of Pa < 10® d, because most of the pro¬ 
genitor systems of short-period binaries produced through 
KL cycles with tidal frictio n initially have Pa < 10® d 
dFabrvckv fc Tremaine! l2007l) . The corresponding semima- 



ioSlo(/shrmk) 


Figure 9. The factor /shrink = “A.i/uA f with which the inner 
binary semimajor axis shrinks for the sampled triples as described 
in Secti 0 n r 4 .ll for isolated triple evolution as in Fig. |8] and here 
binned with respect to the final inner orbital period Paj. Three 
bins in P\ f are assumed, indicated in the legend in the bottom 
panel. The black vertical dotted line indicates the value of /shrink 
chosen for the selection of systems with substantially shrinking 
inner binaries. 


jor axes, a a and ac, are computed from the orbital periods 
using Kepler’s law (neglecting the planet mass m 3 ). 

The eccentricities ba and ec are sampled from a 
Rayleigh distribution, dN/dek oc ej, exp(—/3e^), with rms 
= 0. 33, consistent w ith the results of 
IRaghava n et aL||2010|. In Ragha van et al.] 1I2OIOI) . in contrast 
with IPuauennov fc Mavoil ( 199ll ). the eccentricity distribu¬ 
tion is found not to be substantially different for periods 
Pfe < 10® d and periods Pk > 10® d. Therefore, we choose 
not to sample eccentricities f rom a thermal distribu t ion fo r 
binaries with Pk > 10® d as in iFabrvckv fc Tremainel ll2007l) . 

The inner and outer orbits are assumed to be randomly 
oriented, i.e. their arguments of pericentre ujk and flj, are 
sampled from a random distribution, and their mutual in¬ 
clination lAC is sampled from a distribution that is linear 
in cos(7ac)- Note that, without loss of generality, we fix the 
inner orbit to be aligned with the 2 : axis of the coordinate 
system, i.e. we set ia = 0°. The initial mutual inclinations 
7AB and 7AC are therefore equal to the initial individual in¬ 
clinations Zb and ic, respectively. 
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Figure 10. The distribution of Qo (cf. equation|8ll for the systems 
that exceeded the imposed maximum wall time limit of 15 min. 
The black vertical dotted line shows the value Qo = 1 for which, 
for these systems, planets were assumed to be fully effective (Qo < 
1) or ineffective (Qo > 1) at shielding. 


Sampled triple systems are rejected in the following 
cases at this stage. 

• The stability criterion of iMardling fc AarsethI [ 2 OOII is 
not satisfied. 

• In isolation, the stars in the inner binary would experi¬ 
ence strong tidal interaction, i.e. aA(l — ca) < 3 {Ri + R 2 )- 

• In isolation, the stars in the inner binary would fill their 
Roche lobe, i.e. Rk < Rn.fe for k G {1,2}, where Rh,k is the 
Roche lobe radius computed at p ericentre accord i ng to the 
analytic fitting formulae given bv ISeninskv et ahl |2003). 

2. For the sampled systems of step (1), additional se¬ 
lection is made based on the outcome of the integration of 
isolated systems. This is to ensure that in the absence of the 
planet, the system would remain dynamically stable, and 
that the inner binary is shrunk substantially because of KL 
cycles with tidal friction. The integration time for each sys¬ 
tem is sampled linearly between 1 and 10 Gyr. In addition, 
the following stopping conditions are always imposed. 

• The stars in the inner binary collide or fill their Roche 
lobe, computed as above. 

• One of the three stars evolves past the MS. 


First, the sampled systems are integrated taking into 
account only the inner binary evolution, i.e. stellar and tidal 
evolution. From these systems, we reject those for which the 
stars in the inner binary collide, fill their Roche lobe, or 
evolve past the MS. The former two cases can occur because 
of tidal friction; note that our inner binaries are initially not 
circular, nor synchronized [see Section [5.6.2l and Appendix 
m for discussion on magnetic braking (MB) effects]. 

Secondly, the sampled systems are integrated taking 
into account the evolution of the triple system without a 
circumbinary planet, i.e. taking into account KL cycles in 
the inner and outer binaries, and tidal evolution in the inner 
binary. From these systems, we also reject those for which 
the inner binary stars collide or fill their Roche Lobe. 

In Fig. [7] we show the initial (dotted lines) and fi¬ 
nal (solid lines) inner period distributions for the case 
of an isolated inner binary system (black lines), and 
an isolated triple system (blue lines). Consisten t with 
previous studies on KLCTF (iMazeh fc ShahamI [l97^: 


Eggleto n fc Kiseleva-Eggleton I2OOII : Fabrvckv fc Tremaina 
2OO7I : iNaoz fc Fabrvckvl 12014 1 . the effect of the tertiary 


star is to produce an enhancement of systems with periods 
roughly between 1 and 6 d. 

3. We select systems for which, in the case of an iso¬ 
lated triple, the inner binary shrinks substantially during 
the evolution. To quantify this criterion, we show in Fig. 0 
the distribution of the factor /shrink = CLA,i/cLA,f with which 
the inner binary semimajor axis decreases during the evolu¬ 
tion for the case of isolated triples. Based on this distribu¬ 
tion, we select systems for which /shrink > 1-5. Effectively, 
this implies that the majority of the systems evolve to a sys¬ 
tem with a short-period inner period, 1 < PA,f < 6. This is 
demonstrated in Fig. [51 where we show the distribution of 
/shrink iu the triple case, for different bins of the final inner 
orbital period. 

4. For each of the Ntr = 192 remaining sampled triple 
systems, we define a three-dimensional grid in (ae, 18 ,^ 13 ), 
with size 10 x 10 x 3 = 300. The inclinations ib range from 
0 to 180°, and for the planet mass m 3 we assume three 
values, 10“^, 10“^ and 10° Mj. For ob, we impose the lower 
and upper limits ub,! and hb.u, respectively. 

The lower limit is assumed to be ob,! = 1.5OB.crit.AB, 
where OB.crit.AB is the critical semimajor axis for dynamical 
stability of the planet in a coplanar orbit around the inner bi- 
nary, and which is adopted from the fitting formula given by 
llHolman fc Wiegertlll99f^ . The latter formula is a function 
of the inner binary mass ratio fi = m^lmx and eccentricity 
ca- The factor of 1.5 in ae.i = 1.5OB.crit.AB is a ‘safety fac¬ 
tor’; close to the regi me of dynamical stability, the orbits are 
non-Keplerian (e.g. Lee fc Peald I2OO6I : iLeuing fc Led I2013I : 
iBromlev fc Kenvorj2015l L and, therefore, the orbit-averaged 
approach no longer applies. The upper limit is assumed to be 
ob.u = 0.9aB,crit,BC, where ae.crit.BC is the largest possible 
value of Ob for dynamical stability with respect to the orbit 
of the t ertiary star. The latter val ue is estimated by apply¬ 
ing the IMardling fc Aarseth|l 200 ll criterion to the BC pair, 
with an ‘outer’ mass ratio of gout = m 4 /(mi -|- m 2 -\- m 3 ). 

In the analysis below, we use the dimensionless quantity 
as which is closely related to ob, and is defined as 


CtB = 


HB — CLB,l 
flB.u — flB.l 


( 10 ) 
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Figure 11. The distributions of the parameters as (first column), Qo (second column), ib (third column) and m 3 (fourth column) for 
systems in which the orbit of the planet remains stable during the entire evolution, and no other stopping conditions occurred. Systems 
are binned with respect to the final inner orbital period each row corresponds to a different bin, indicated in the first column. Also 

indicated in the first column is the number of systems in the corresponding bin of 


This quantity is motivated by the fact that the allowed range 
of initial as for dynamical stability varies per triple system. 
By construction, this is not the case for qb, which, for any 
sampled triple, ranges between qb = 0 (ub = 1.5UB.crit.AB, 
i.e. the planet’s orbit is close to the dynamical stabil¬ 
ity limit with respect to the inner binary), and ob = 1 
(ub = 0.9 UB.crit.BC, i-s. the planet’s orbit is close to the dy¬ 
namical stability limit with respect to the outer binary). In 
other words, by using an, we can evaluate the effect of the 
circumbinary planet for a population of triples with different 
parameters. 

For the A^quad = 300 x 192 = 57600 integrations with a 
planet, the following additional stopping conditions are set. 

• The circumbinary orbit is unstable with respect 


to the inner and oute r binaries according to the 
llMardling fc Aarsethll200l] ) criterion. 

• The circumbinary orbit intersects with the inner binary, 
i.e. aB(l — cb) < oa(1 + ba), and is therefore unlikely to be 
stable. 

For some quadruple systems, the integration proceeded 
slowly. We therefore imposed a maximum wall time of 15 
min per system, which was reached for 19% of the inte¬ 
grated systems. In Fig. 1101 we show the distribution of Qo 
(cf. equation [S]) for the systems in which the maximum wall 
time was exceeded. For ~ 60% of these systems, Qo < 1. 
As shown previously in Section[3l for Qo < 1, shielding is 
typically effective, whereas it is not for Qo > 1. Therefore, 
for the purposes of the analysis below, we will assume that 
in the systems that exceeded the integration wall time, if 
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2 o 1 , shielding is completely effective and hence the in¬ 
ner orbital period does not change (i.e. PA,f = PA,i), and 
if Qo ^ 1 , shielding is completely ineffective and hence the 
inner orbital period shrinks as much as would be the case 
without a planet. Clearly, this is approximate because in re¬ 
ality, there is no sharp transition between the two regimes 
and there is also a dependence on the inclination ie (cf. 
Fig.iJ. However, we do not expect that the approximation 
strongly affects our conclusions. 


4.2 Results 

4-2.1 Initial parameters for surviving planets as a 
function of the final binary period 

In Fig. 111! we show the distributions of the initial parame¬ 
ters Ob (cf. equation [Toll . Qo (cf. equation [ 8 |, is and m 3 in 
the first through fourth columns, respectively, for systems 
in which the orbit of the planet remains stable during 
the entire evolution according to the criteria discussed in 
Section im Each row in Fig. [11] corresponds to a different 
bin of the final inner orbital period, PA,f- 

For short periods, the planet tends to be further away 
from the inner binary (i.e. qb and Qo are typically larger) 
compared to longer periods. This trend is compatible with 
the lack of tight coplanar planets around short-period bina¬ 
ries, although for our sampled population, there is certainly 
no absolute lack of planets in tight (i.e. qb < 0.3) orbits 
around short-period binaries. 

In fact, for Id < Pa,! < 6 d, the distribution of as 
is still peaked around as ~ 0.1. There are two important 
factors that contribute to this. 

(i) In our initial triple population, there is a non- 
negligible number of systems with initial inner periods 1 d < 
Pa,! < 6 d (cf. Fig. [7|. Evidently, for these systems, a pos¬ 
sible outcome is a final period 1 d < Pa,! < 6 d when the 
inner orbit did not shrink due to effective shielding by the 
planet, likely corresponding to a massive planet in a tight 
and/or coplanar orbit. In other words, in principle, a short- 
period binary can simply be born with a tight and coplanar 
planet, and the planet would prevent the inner orbital pe¬ 
riod from becoming even shorter. Given the observational 
evidence for a t hird stellar companion around most short- 
period binaries jTokovinin et al.ll2006ll . one might consider 
the possibility that such binaries can not primordially form 
without the effects of KLCTF. If this is the case, then the 
primordially short-period binaries in our simulations should 
not be considered, and therefore, we would expect no tight 
coplanar circumbinary planets to exist at all around short- 
period binaries. 

(ii) For Ob // 0.3 or Qo ^ 10°, the planetary orbit for a 
large fraction of systems becomes highly eccentric and inter¬ 
sects with the inner binary, and likely becomes unstable (cf. 
Sections l4.2.2l and 15.31) . These cases are not included in Fig. 
El and this produces a bias for an absence of stable systems 
with Qb ^ 0.3. 

To take these complicating factors into account, we in¬ 
troduce the ‘shielding efficiency’ ? 7 shieid which quantifies the 
planet’s ability to shield the inner binary from KL eccen¬ 


tricity oscillations induced by the tertiary star, defined as 

_ /"shrink,triple /"shrink,quad i \ 

^shield = 7 Z • 

/shrink,triple J- 

Here, /shrink.tripie is the factor with which the inner binary 
semimajor axis decreases in the absence of a planet, and 
/shrink,quad is the Corresponding factor with a planet present. 
Note that we restrict to systems for which /shrink.tripie > 
1.5. If the planet is able to fully shield the inner binary, 
/shrink,quad = 1, aud T^shieid = 1- Ou the Other hand, when 
shielding is completely ineffective, /ahrink,quad = /shrink.tripie, 
and J^shieid = 0. For some systems, we find that the planet 
can enhance the inner binary eccentricy excitations; in this 
case, ?7shieid < 0 (cf. Section l4.2.4l) . 

In Fig. 1121 we show the distributions of ob, Qo, *b 
and m 3 , where the systems are now binned with respect 
to Tjshieid. "When shielding is ineffective (r^shieid ~ 0, cf. the 
second row in Fig. El), Qo is invariably > 1 and m 3 tends 
to be low (i.e., the large value of m 3 = 1 Mj is disfavoured). 
In contrast, when shielding is highly effective (r^shieid ~ 1, 
cf. the last row in Fig. El), Qo is typically < 1 (with a peak 
near Qo = 0.1), and m 3 tends to be high (i.e., the large 
value of m 3 = 1 Mj is favoured). 

The distribution of the initial in in Fig. [12] is a strong 
function of ob and m 3 . In Fig. 1131 we show, for 0.9 < 
JJshieid < 1 , the inclination distribution for different ranges 
of Ob and m 3 . For ob < 0.3, there is a tendency for coplanar 
orbits, for both low-mass and massive planets. For planets 
in wider orbits, ob > 0.3, the distribution of in is strongly 
peaked towards ib ~ 90° for low planet masses, whereas 
for high planet masses, the distribution of in is much less 
peaked. These trends are qualitatively consistent with the 
smaller number of integrations that were carried out in Sec- 
tion[ 3 ]for fixed parameters of the stellar triple. 

In Fig. [14] we show, for each of the final period bins 
of Fig. 111 ! the initial distributions of qb, Qo, *b and m 3 , 
where additional binning was made with respect to » 7 shieid- 
Black (blue) lines correspond to ryshieid < 0.2 (? 7 shieid > 0.2). 
Cumulative distributions for the different final period bins 
and the two ranges of shielding efficiencies are shown in Fig. 
E] As expected, the close planetary orbits from the top row 
in Fig. E] (Id < Pa,! < 6 d) typically correspond to a high 
shielding efficiency (cf. the blue lines in Fig. I14II . This indi¬ 
cates that the inner binary did not shrink, but was simply 
formed with a planet in such an orbit. On the other hand, 
the wider orbits typically correspond to a low shielding effi¬ 
ciency (cf. the black lines in Fig. I14II . In the latter case, the 
initial inner orbit was much wider, and shrank substantially. 

4-2.2 Unstable planets 

As mentioned in Section l4.2.1l for a large fraction of sys¬ 
tems, the planetary orbit becomes highly eccentric because 
of the torque of the outer orbit, and intersects with the in¬ 
ner binary. This likely results in an unstable planetary orbit, 
with a high probability for the planet being ejected from the 
system (for further discussion on and calculations of the out¬ 
comes, see Section lS^ . Such orbit crossings typically occur 
early in the evolution. This is illustrated in Fig. 1161 where 
we plot the distribution of the time of orbit crossing Gross, 
when applicable. For ~ 80% of the cases, tcross < 10 Myr, 
which is very short compared to the MS time-scale of any of 
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Figure 12. The distributions of the parameters ckb (first column), Qo (second column), zb (third column) and m 3 (fourth column), 
where systems are binned with respect to the shielding efficiency Shield (cf- equation IIID : each row corresponds to a different bin, 
indicated in the first column. Also indicated in the first column is the number of systems in the corresponding bin of z?shield- 
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Figure 14. The distributions of the parameters aa (first column), Qo (second column), ib (third column) and m 3 (fourth column) for 
systems in which the orbit of the planet remains stable during the entire evolution, and no other stopping conditions occurred. Systems 
are binned with respect to the final inner orbital period Pji^ f as in Fig. 1111 here, we also make a distinction between a low shielding 
efficiency () 7 shield < 0 - 2 ; black lines) and high shielding efficiency (r^ghield > 0 - 2 ; blue lines). 


the stars (as described in Section ldTIl the integration time 
is sampled linearly between 1 and 10 Gyr). Because of this 
quick removal of the planet, it is unlikely that the planet 
could affect the inner binary evolution in these cases, and, 
therefore, we correspondingly assume tishieid = 0. 

In Fig. 1171 we show the fraction /cross of systems with 
orbit crossing as a function of the initial parameters qb, Qo, 
iB and m 3 . This fraction increases strongly from /cross ~ 0.4 
for ob ~ 0 to /cross ~ 0.9 for as > 0.5. This may seem coun¬ 
terintuitive: for larger ob, the planet is further away from 
the inner binary, therefore orbit crossings are expected to be 
less likely. However, for larger qb, the planet is also placed 
closer to the tertiary star, causing greater eccentricity ex¬ 


citation in the circumbinary orbit (e.g. the top-left panel 
in Fig. [5|. Apparently, the latter effect is (much) stronger 
than the former, resulting in more orbit crossings for wider 
planetary orbits. As mentioned in Section [4.2.1l the high oc¬ 
currence of orbit crossings produces a tendency for a lack of 
stable planets for ob ^ 0.3. 

There is no strong dependence of /cross on is, nor m 3 . 
This can be understood by noting that the distribution of 
the mutual inclination between the circumbinary and the 
outer orbits, isc, is the same as the distribution of the 
mutual inclination between the inner and outer orbits, iac 
( i.e. cos[iAc] was sampled linearly between -1 and 1; cf. Sec- 
tion l4.ll) . Therefore, the distribution of isc is independent 
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Figure 13. The distributions of the initial ib for 0.9 < )7shield < 1 
(cf. Fig. I12II . distinguishing between different ranges of «b a-nd 
m 3 . In the top (bottom) row, ae < 0.3 (ae > 0.3); in the left 
(right) column, m 3 = 10“^Mj (m 3 = lO'^Mj). The number of 
systems in each case is indicated in each panel. 


of iB, and there is no preference with regard to /cross for 
orbits with ib close to 90°. Furthermore, because the planet 
mass m 3 is typically less than 10 “® of the stellar masses, 
Fkl.bc oc (nil + m2 + m 3 + m4,)lm4, is nearly independent 
of m 3 , explaining why /cross is essentially independent of 
m 3 . Note that, in contrast, the shielding efficiency tjshieid is 
a strong function of m 3 (cf. Fig. 1121) . 


4-2.3 Final planet orbital eccentricities 

In Figs 1181 and 1191 we show the distributions of the planet 
orbital eccentricity at the end of the integration, CB.f, binned 
with respect to PA,f and r^shieid, respectively. In the case of 
binning with respect to T^shieid, orbit crossing cases are also 
included, in which case we set r^ahieid = 0 (cf. Section l4.2.2ll . 
For a crossing of the planetary orbit with the inner binary to 
occur, cb needs to be high. This causes CB.f to be typically 
high if ryshieid « 0; in « 50% of the cases, en.f > 0.9. On the 
other hand, if the planetary orbit remains stable, cb tends 
to be much smaller. For cases of very effective shielding, i.e. 
0.9 < ryshieid < 1, ~ 90% of the systems have en.f < 0.2. In 
Fig. 1181 the planetary orbit remains stable in all cases, and 
eB,f tends to be low. There is little to no dependence of the 
distribution of cbj on the final orbital period. 


4 . 2.4 Reversed shielding 

As shown in the top row of Fig. 1121 in a relatively small 
number of cases, r^ahieid < 0 , i.e. the inner binary shrinks 
more compared to the situation without a planet. The val¬ 
ues of ctb for which this occurs are strongly peaked towards 
small values, and Qq is peaked near Qo = 1- For Qo ~ 1, the 
KL time-scales for the AB and AC pairs are approximately 
equal, and, as similarly mentioned with respect to the quan¬ 
tity TZo by HPAPZ15, this can give rise to chaotic dynamics 
and possible enhancement of the eccentricity in the inner bi¬ 
nary, and hence a negative r?shieid. Although interesting from 



Figure 15. Cumulative distributions for Qo, binned with re¬ 
spect to the final inner orbital period Pa f (indicated in the leg¬ 
end in the top panel). In the middle (bottom) panels, a distinc¬ 
tion is made with respect to the shielding efficiency r^gjjieiu (cf. 
equation Hil l: including the systems with »?shield < 0.2 (middle 
panel), and with r)sj,ieiu > 0.2 (bottom panel). The meaning of 
the linestyles is the same between the different rows. 


a theoretical dynamical point of view, negative ? 7 shieid only 
occur for a small number of systems, i.e. a fraction of 0.032 
of the stable systems, and 0.0084 of all systems. 
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Iogio(m3/Mj) 


Figure 17. The fractions of systems for which the planet either remains stable (blue solid lines), or for which the orbit of the planet 
becomes highly eccentric and intersects with the inner binary, likely resulting in an unstable orbit (red dashed lines), as a function of 
So, *B and m3. 


5 DISCUSSION 

5.1 An approximate analytic condition for planet 
shielding — implications for other systems 

In Section ISTi we have shown that planet shielding is effec¬ 
tive for Qo ^ 1 (cf. Fig. [2]). We have also confirmed this 


for a much larger number of systems, and with the inclusion 
of tidal friction, in Section|4] (cf. Fig. I12II . As demonstrated 
in e.g. Figs [2] |4] and m there is also a dependence of the 
planet shielding ability on the inclination ib. Nevertheless, 
as a first approximation, the condition Qo ^ 1 can be used to 
evaluate the importance of planet-shielding in stellar triples. 
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logio(^ cross /Myr) 


Figure 16. The distribution of the time of orbit crossing tcrosa 
of the planet with the inner binary, for the systems in which this 
occurred. 


Using equation (( 8 |, the condition can be written as 
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1 “ ec,o 
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1 "^B.D 
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where in the second line, we assumed that m 3 is negligible 
compared to the other masses, and in the last line, we as¬ 
sumed cb.o = ec,o = 0. For small planet masses, the scaling 
with m 3 is Ob oc mj^^. Equation [T^] implies that for a planet 
of order Jupiter mass in a stellar triple (i.e. ruk ~ 1 M© 
for k G {1,2,4}), the circumbinary semimajor axis should 
be less than approximately a tenth of the outer semimajor 
axis. 

We remark that the fourth body orbiting the circumbi¬ 
nary planet system may also be less massive (by factors up 
to a thousand) than was assumed in Sections[3] and [d] In 
particular, it could be a brown dwarf or a massive planet 
orbiting a lower-mass circumbinary planet. As decreases, 
equation m implies that for fixed ac, ub can be larger for 
shielding to be effective. This can be explained intuitively 



Figure 18. The distributions of the final eccentricity 63 ,f of the 
planetary orbit for stable systems and binned with respect to the 
final inner orbital period P\ f. 


by noting that for fixed m 3 and decreasing m 4 , the Pkl,ab 
time-scale remains constant, whereas the Pkl.ac time-scale 
increases, i.e. the outermost body becomes less dominant. 
This, combined with the lower efficiency of a tertiary with a 
lower mass to shrink the inner orbit, suggests that the lack 
of circumbinary planets around short-period binaries is even 
more severe for triples with low-mass tertiary companions. 
Similarly, shielding is expected to be more effective for lower 
inner binary masses mi and m 2 , suggesting a more severe 
lack of planets around low-mass short-period binaries. 


5.2 Shielding of the planetary orbit by the inner 
binary 

Although not discussed in detail here, shielding can also oc¬ 
cur in the orbit of the planet: the inner binary can induce 
rapid precession in the planetary orbit, shielding the latter 
against high-amplitude KL eccentricity oscillations induced 
by the outer binary. This aspect was discussed in detail by 
HPAPZ15, where it was shown that the approximate con¬ 
dition for shielding of the planetary orbit is TZo ^ 1, or (cf. 
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become highly eccentric because of KL eccentricity oscilla¬ 
tions induced by the tertiary companion, and intersect with 
the inner binary. In fact, this is a likely scenario, as demon¬ 
strated by e.g. Fig. 1171 Here, we investigate the possible 
outcomes of such orbit crossings and hence likely unstable 
orbits, by integrating the fo ur-body system using the di¬ 
rect A'^-body code Hermite jHut et al. 1995ll incorporated 
in A MUSE dPorteeies Zwart et al.l ^13l : Pelupessv et al.l 
I 2 OI 3 II . Here, tidal effects are not included. The initial con¬ 
ditions are taken from the grid, of size 300, of the initial 
parameters ob, *b and m 3 for one the triple systems sam¬ 
pled in Section|31 Details of the parameters are given in the 
third column of Table [T] 

First, the systems are integrated using our secular code 
until the circumbinary orbit intersects with the inner binary. 
Subsequently, we sample 100 different sets of random mean 
anomalies for the three orbits, and integrate the system for 
each set for the duration of 40 initial circumbinary orbital 
periods. In this manner, the most important outcomes are 
revealed for each combination of planet parameters. In total, 
300 X 100 = 3 X 10"* direct A'^-body integrations were carried 
out. 

We find the following outcomes, in order of decreasing 
likelihood /: 

(i) the planet becomes unbound from the system 
if ~ 0.657); 


Figure 19. The distributions of the final eccentricity eB,f of the 
planetary orbit, binned with respect to the shielding efficiency 
Shield- Here, cases are included when the planetary orbit inter¬ 
sects with the inner binary, in which case r?shield is assumed to be 
0 . 


equation [9]) 


^ 1/3 2/3 

^ CL A CLq 

\m4j 
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where in the second and last lines, we neglected m 3 com¬ 
pared to the other masses, and where we assumed circu¬ 
lar orbits. In general, TZq ';$> Qo (cf. equation [9]), indicating 
that typically, shielding of the inner binary by the planet 
(Qo ^ 1) is more likely than shielding of the planet by the 
binary (TZo < 1). 

We note that similar dynamics apply to satellites 
around the Pluto-Charon binary sy stem in the Solar system. 
As shown bv iMichaelv et akl (l2015l l. if the orbit of a satellite 
is close enough to the Pluto-Charon binary, then precession 
induced by the Pluto-Charon binary on the satellite’s orbit 
protects the latter from KL oscillations induced by secular 
perturbations of the Sun. 


5.3 The fate of planets with unstable orbits 

As mentioned in Sections[3] and [4l when the orbit of the 
planet is relatively close to the outer binary, the former can 


(ii) the planet orbits the outer binary (i.e. a circumtriple 
planet; / « 0.289) 

(hi) the planet collides with a star (/ ~ 0.029); 

(iv) the planet remains stable as a circumbinary planet, 
but with a different orbit (/ « 0.017); 

(v) the planet becomes bound to a single star (/ ~ 0.009). 

The large probability of ejection of the planet is intuitively 
easy to understand from the low mass of the planet com¬ 
pared to that of the stars. Circumtriple planets are also com¬ 
mon, but note that the integration time is limited and the 
orbits are typically eccentric, therefore not all systems may 
be stable indefinitely. 

In the top row of Fig. 1201 we show the distribution of 
the fractions of the outcomes. In the bottom row of the same 
figure, we show the distributions for collisions with stars 
(left column) and bindings to stars (right column), where 
a distinction is made between the different stars. Collisions 
are most common with the primary star. This is easily un¬ 
derstood by noting that by definition, the primary star is 
the most massive, and, therefore, also has the largest radius 
(as described in Section[^ the radii are calculated using the 
SeBa stellar evolution code). With regards to the planet 
becoming bound to a single star, being bound to the ter¬ 
tiary is most likely. This may be the result of a lower orbital 
speed of the planet when it is close to the tertiary (roughly 
corresponding to the apocentre of the circumbinary orbit if 
it were still stable), as opposed to when it is close to the 
inner binary (roughly corresponding to the pericentre of the 
circumbinary orbit if it were still stable). 

The outcomes found above - keeping in mind the caveat 
that they are based on integrations of only a single triple sys- 
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tern and that the planetary parameters were taken from a 
linear grid - have interesting implications for stellar triples. 
Case (i), combined with the high fraction of orbit-crossing 
planets as found in Section [4.2.2l suggests that circumbi- 
nary planets in triples are likely to become unbound from 
their parent binary early in their evolution, and become free- 
floating planets. Case (ii) suggests that circumtriple planets 
should be fairly common. Case (hi) provides a scenario for 
polluting stars in the inner binary, the primary star in par¬ 
ticular. Lastly, case (v) suggests that in stellar triples, cir- 
cumbmary planets (i.e. P-type planets) can be transformed 
into circumste//ar planets (i.e. S-type planets), albeit with a 
low probability. 


5.4 Implications for planets around blue straggler 
stars 


As first suggested by IPerets fc Fabrvckvl (l20 09h. a n d late r 
studied more quantitatively by Naoz fc Fabrvck'^ (l2014ll . 
KLCTF may also lead to mergers or mass transfer in short- 
period binaries, producing blue straggler stars (BSSs). Our 
results suggest (secular dynamical) constraints for planets 
around BSSs formed in this manner, where the planets are 
either in a circumbinary configuration (in case of binary 
BSSs) or in a circumstellar configuration (in case when the 
original stellar binary merged to produce a single BSS star). 
Analogously to the case of short-period binaries, we expect 
such planets around BSSs to be typically of low mass, and in 
wide and/or inclined orbits around the BSS or BSS binary. 


5.5 Approximations in the integrations 

In the numerical integrations of Sections[3] and [4] the ‘cross’ 
term that appears in the Hamiltonian at the octupole-order, 
Hoct,cross (cf. section 2.4 of HPAPZ15), was neglected. Note 
that the ‘non-cross terms’, Lfoct.AB, Hoct,BC and Hoct.AC, 
were always included in Sections[3] and ID For the systems 
considered here, it is unlikely that the cross term Lfoct,cross 
has a large effect on the dynamical evolution because its nu¬ 
merical value is generally very small compared to the other 
terms that appear in the octupole and the next higher, hex- 
adecupole, orders. For example, for the system chosen in 
Section lSTT] (cf. Table [1] setting ms = 1 Mj, f b = 0 ° and 
ic = 90'’), the ratio r of the absolute value of the orbit- 
averaged cross term to the absolute value of all other orbit- 
averaged terms at octupole and hexadecupole order, dehned 
in equation (10) of HPAPZ15, is r ~ 6 x 10“® if ub = 1AU, 
and r « 7 X 10“^ if ob = 5 AU. For the system of Section[3]2] 
(setting ms = 1 Mj, is = O'’ and ic = 90'’), r « 6 x 10“^^ if 
Ob = 10 AU, and r « 2 x 10~^ if ob = 50 AU. 

Furthermore, in the integrations we assumed that there 
is only one planet orbiting the inner binary. Among the 
currently limited number of Kepler transiting circumbinary 
planets, there is already one system, Kepler 47 (cf. Table 
Hr]) with three confirmed circumbinary planets orbiting a 
single binary. Although the case of multiple planets is be¬ 
yond the scope of this work and left for future work, we note 
that in this situation, the planet with the tightest orbit with 
respect to the inner binary likely has the largest potential 
for shielding the inner binary from the secular torque of the 
outer orbit. This picture is complicated by the possibility of 


planet-planet scattering and mean motion resonances if the 
planets are in close orbits to each other. 


5.6 Other dissipative effects 

In this work, we considered the dissipative effect of tidal fric¬ 
tion in the inner binary. Other potentially important dissi¬ 
pative processes that affect the orbital energies are gravita¬ 
tional wave (GW) emission, MB and gas drag in circumbi¬ 
nary discs. 


5.6.1 Gravitational wave emission 

Shrinkage due to GW emission was not included for nu¬ 
merical reasons, but we note that it is important only in 
very tight binaries that are beyond the scope of this pa¬ 
per. For circular binaries, the GW inspiral time is tew = 
Sc^Ua /f256G’ ^mim2(mi -I- m 2 )], where c is the speed of light 
(|Peterslll964l L which is equal to 10 Gyr if logj^Q(PA/d) « 
—0.4. In our population synthesis, the final orbital periods 
are typically longer than this value (cf. Fig. [7]), indicating 
that shrinkage in the inner orbit due to GW emission is 
not important (note that the inner orbit is typically nearly 
circular when the final orbital period is reached). 


5.6.8 Magnetic braking 

Many uncertainties still exist regarding the efficiency of MB. 
We evaluated the effect of MB by carrying out the integra¬ 
tions of isolated binaries and triples as in Section lUTI now 
also including MB for the primary and secondary stars in 
the inner binary, with standard assumptions for the MB law 
and its efficiency. Further details are given in Appendix IbI 
We find that when including the effects of MB, there is 
a nearly absolute absence of binaries with periods < 10 d 
(cf. Fig. \A^, which is grossly inconsistent with observa¬ 
tions. The latter, in contrast, show a peak in the distribu¬ 
tion around 3-6 d llTokovininll2014h . This suggests that MB 
is likely not as efficient as was assumed here. Because MB 
is not the focus of this work, we chose to exclude it in the 
integrations of Sections[3] and |4l Nevertheless, the apparent 
discrepancy between the canonical efficiency of MB and the 
period distribution of solar-type MS binaries merits future 
investigation. 


5.6.3 Drag in circumbinary gas discs 

If there is a disc around the inner binary, gas drag can cause 
the circumbinary planet to migrate inwards. Consequently, 
Qo would decrease (cf. equation [SJ , and hence the planet’s 
shielding ability would be increased. Howe ver, the lifetime o f 
circumbinary discs is at most a few Myr llAlexanderj|2012| L 
which is much shorter than the typical duration of our inte¬ 
grations (in our population synthesis, the integration time 
varied between 1 and 10 Gyr to reflect the typical age of ob¬ 
served Kepler-like systems). Therefore, disc migration likely 
does not affect our conclusions regarding planet shielding in 
Kepler-like systems. 
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Figure 20. The distributions of the fraction of direct Af-body in¬ 
tegrations with various outcomes, in the case that the orbit of the 
planet becomes highly eccentric and intersects with the inner bi¬ 
nary. The triple system is fixed, whereas , *b and m 3 are taken 
from a grid with size 300 as in Section[4!T] Top row: overview, 
with the major channels described in the legend. Bottom row: 
making a distinction between stars for cases when the planet col¬ 
lides with a star (left column), or when it becomes bound to a 
star (right column). 


6 CONCLUSIONS 

In recent observations of circumbinary planets, there is an 
unexpected lack of coplanar circumbinary planets around 
short-period solar-type MS binaries. This goes against ob¬ 
servational biases of detecting such planets. We have shown 
that the lack can be explained by the secular gravitational 
influence of a circumbinary planet in hierarchical triple sys¬ 
tems. Our arguments and conclusions are as follows. 

1. Observations show that binaries with periods of ~ 3 — 6 d 
are nearly always (96% ) orbited by a tertiary companion star 
dTokovinin et al.ll2006l ). The short orbital period likely re¬ 
sulted from Kozai-Lidov (KL) eccentricity oscillations in the 


inne r orbit induced by the tertiary, combined with tidal fric- 


tion (lMazeh_&_Slmha^[l97^ ^deton ^JUselcvn^Eeglet^ 

I 2 OOII : IPabrvckv fc Trernaindl200^Naoz fc Fabrvckvl 20141 ). 


This suggests that many of the short-period Kepler eclipsing 
binaries are triple star systems. The progenitor inner binary 
was likely wider, with an orbital period of up to ~ 10 “* d. 


2. We have demonstrated that if there is a circumbinary 
planet around the progenitor inner binary, then the KL ec¬ 
centricity oscillations in the inner orbit induced by the ter¬ 
tiary can be quenched through the secular gravitational in¬ 
fluence of the circumbinary planet. Thereby, the inner binary 
is ‘shielded’ from the torque of the tertiary star, and does 
not shrink to a tight orbit. However, this only occurs if the 
circumbinary planet is sufficiently massive, and if its orbit 
is sufficiently close to and coplanar with the inner binary. In 
many other cases, the circumbinary orbit is stable but can¬ 
not efficiently shield the inner binary, or the circumbinary 
orbit is destabilized and is most likely ejected. 


3. In particular, if a low-mass circumbinary planet is initially 
inclined with respect to and far away from the inner binary, 
then planet shielding is typically ineffective, and the inner 
binary can shrink to a tight orbit. On the other hand, for a 
more massive planet in an initially approximately coplanar 
and tight circumbinary orbit, shielding typically prevents 
the inner orbit from shrinking (cf. Figs ll4l and ll2ll . 

4. Consequently, for systems with circumbinary planets sur¬ 
viving through the evolution, short-period inner binaries 
typically do not have massive (m 3 ~ Mj) circumbinary plan¬ 
ets on tight (Qo < 1 , where Qo is defined in equation[Sl) and 
coplanar orbits. Namely, if the latter was present initially, 
the inner binary would unlikely have shrunk to a short pe¬ 
riod. In contrast, if the circumbinary planet is less massive 
(m 3 ~ 10~^ Mj) and initially inclined with respect to and 
far away from the inner binary (Qo > 1), then shrinkage 
is possible. This trend is consistent with the current Kepler 
observations. The transition with respect to the planet mass 
and orbital semimajor axis occurs at Qo ~ 1, which, for a 
Jupiter-mass planet in a solar-mass triple, corresponds to a 
circumbinary semimajor axis roughly equal to a tenth of the 
outer orbit semimajor axis (cf. eauation ll2D . 

Our results suggest that similar constraints also apply 
to planets around blue straggler stars. 
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Figure Al. The initial (black) and final (blue) inner period dis¬ 
tributions for the sampled triples as described in Section r4.1l (cf. 
Fig. ETJ, here also including the effects of MB as discussed in 
Section r5.6.2l and Appendix[3 Distributions corresponding to iso¬ 
lated binary (triple) evolution are shown with black (blue) lines. 
The initial distributions are shown with dotted lines. Cases when 
the inner binary merges are shown with dashed lines. Otherwise, 
the final distributions are shown with solid lines. 


APPENDIX A: KEPLER TRANSITING 
CIRCUMBINARY PLANETS 

In Table HH we give an overview of the currently known 
Kepler transiting circumbinary planets. 



0.0 0.5 1.0 1.5 2.0 

loglO (/shrink) 


Figure A2. The factor /shrink = aA,i/iiA,f with which the in¬ 
ner binary period shrinks for the sampled triples as described 
in Secti 0 n r 4 .ll here also including the effects of MB as discussed 
in Section r5.6.2l and Appendix]^ Only systems are included for 
which the inner binary does not merge. Black (blue) lines: iso¬ 
lated binary (triple) evolution. The black vertical line indicates 
the cutoff value chosen for our population synthesis study (cf. 
Section l4.1ll . 


APPENDIX B: MAGNETIC BRAKING IN 
TRIPLES 

As discussed in Section r5.6.2l we carried out the integra¬ 
tions of isolated binaries and triples as in Section lTTI now 
also including the effects of MB for the primary and sec- 
o ndary stars in the inner binary by adopting equation (1) 
of I Barker fc Ogilvid ll2009l). The latter eq uation is based on 
the Skumanich rel a tion |Skumanicblll972h and the re sults of 
IVerbunt fc ZwaanI lll98ll) : IPobbs-Dixon et all (|2004| ). Here, 
we assumed a MB pa rameter of omb = 1.5 x 10”^^ yr 
llBarker fc Ogilvijl2009h . 

In Fig. we show the distribution of the initial and 
final inner orbital periods as in Fig. [7] but now with MB 
included. MB has a large effect on the orbital evolution. In 
contrast to the situation without MB, many systems merge 
during the evolution even for isolated binary evolution (cf. 
the dashed lines in Fig. mj. Furthermore, when including 
the tertiary star, no peak is produced in the dual period 
distribution for Pa,! ~ 1 — 6d as in Fig. [7] Instead, owing 
to the large number of merged systems, there is an absence 
of systems with these periods. 

For the non-merging systems, we show the distributions 
of semimajor axis shrink factor /shrink in Fig. IA2I Although 
the inner binary still shrinks more in the isolated triple case 
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System 

mi 

m2 

m3 

Ri 

R2 

Rs 

a-A 

Pa 

EA 

as 

Pa 

eb 

Reference 


M© 

Mq 

Mj 

R© 

R© 

Rj 

AU 

d 

- 

AU 

d 

- 


Kepler 16b 

0.690 

0.203 

0.333 

0.649 

0.226 

0.754 

0.224 

41.1 

0.159 

0.705 

229 

0.0069 

a 

Kepler 34b 

1.05 

1.02 

0.220 

1.16 

1.09 

0.764 

0.229 

27.8 

0.521 

1.09 

289 

0.182 

b 

Kepler 35b 

0.888 

0.809 

0.127 

1.03 

0.786 

0.728 

0.176 

20.7 

0.142 

0.603 

131 

0.042 

c 

Kepler 38b 

0.949 

0.249 

<0.384 

1.76 

0.272 

0.388 

0.145 

18.8 

0.103 

0.464 

106 

<0.032 

d 

Kepler 47b 

1.04 

0.362 


0.964 

0.351 

0.266 

0.0836 

7.45 

0.0234 

0.296 

49.5 

<0.035 

e 

Kepler 47c 

1.04 

0.362 


0.964 

0.351 

0.411 

0.0836 

7.45 

0.0234 

0.989 

303 

<0.411 

e 

Kepler 47d 

1.04 

0.362 


0.964 

0.351 


0.0836 

7.45 

0.0234 

0.72 

187.3 


f 

Kepler 64b 

1.47 

0.37 


1.7 

0.34 

0.55 

0.177 

20.0 

0.204 

0.642 

139 

0.1 

g 

Kepler 413b 

0.820 

0.542 

0.211 

0.776 

0.484 

0.387 

0.101 

10.1 

0.0365 

0.355 

66.3 

0.118 

h 

KIC 963289 

0.934 

0.194 

3 X 10 -* 

0.833 

0.214 

0.549 

0.185 

27.3 

0.0510 

0.788 

241 

0.0379 

i 

ROXs 42B b 

0.89 

0.36 

6-14 







~140 



j>k 

FW Tau b 

0.28 

0.28 

6-14 







~330 



j,k 

Ross 458 C 

0.6 

0.09 

6-12 







150-800 



k 

SR12 C 

1.05 

0.5 

12-15 







~ 1100 



k,l 

DP Leonis 

0.6 

0.09 

6.39 




0.027 

0.062 

0 

8.2 

28 

0.39 

m 


Do^le^^aU 


Welsh et al, 


Welsh et al 


Orosz et al 



Qrosz et al. 


^ Orosz (in prep) 


s Rostov et al. 

20131 

^ Rostov et al. 

(2014) 

^ Welsh et al.| 1 ^ 

ioisll 

ioill 

^ Kraus et al.l 

^iBowler et al. 

O 

S- 

^Kuzuhara et al.l (1201 iT) 

"^ISchwooe et a 

1.1 ll2002h:l0 


ian et al 


euermann 


et al. Hm3) 


Table Al. Observed circumbinary planets. The binary primary and secondary masses (radii), mi (-Ri) and m 2 (R2), respectively, are expressed 
in Mq (R 0); the circumbinary planet mass (radius), m3 (-R3), is expressed in Mj (Rj). The binary and circumbinary orbits are denoted with A 
and B, respectively. The semimajor axes, a\ and as, are expressed in AU, and the orbital periods. Pa and Pb, are expressed in days. Note that 
so far, three circumbinary planets have been detected around the binary in the Kepler 47 system. Kepler 64b is also known as KIC 4862625; the 
binary+planet system is orbited by another stellar binary. Top rows: Kepler transiting circumbinary planets. Bottom rows: circumbinary planets 
detected through other methods. 


compared to the isolated binary case, the enhancement in 
/shrink IS not nearly as large compared to the situation with¬ 
out MB (cf. Fig. [8]). In other words, MB is so efficient at 
driving an inner binary merger in isolation, that KL cycles 
with tidal friction hardly affect the evolution. 

As also discussed in Section [5.6.2l this suggests that MB 
is likely not as efficient as was assumed here, and this ap¬ 
parent discrepancy between the canonical efficiency of MB 
and the period distribution of solar-type MS binaries merits 
future investigation. 
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